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^ I Abstract 

. We address some current theoretical issues around ultra-high energy cosmic rays. We recall 

' that scenarios producing more 7— rays than cosmic rays up to high redshift can in general only 

. provide a sub-dominant contribution to the ultra-high energy cosmic ray flux. This includes 

' extra-galactic top-down and the Z-burst scenarios. Finally we discuss the influence of large 

. scale cosmic magnetic fields on ultra-high energy cosmic ray propagation which is currently 

' hard to quantify. The views presented here represent the authors perspective. 

P ; 1 Introduction 

! High energy cosmic ray (CR) particles are shielded by Earth's atmosphere and reveal their existence 

^ ' on the ground only by indirect effects such as ionization and showers of secondary charged particles 

■ covering areas up to many km^ for the highest energy particles. In fact, in 1912 Victor Hess 

, discovered CRs by measuring ionization from a balloon [T], and in 1938 Pierre Auger proved the 

' existence of extensive air showers (EAS) caused by primary particles with energies above 10^^ eV 

by simultaneously observing the arrival of secondary particles in Geiger counters many meters 

apart 0. 

After almost 90 years of research, the origin of cosmic rays is still an open question, with a degree 
of uncertainty increasing with energy Only below 100 MeV kinetic energy, where the solar wind 
shields protons coming from outside the solar system, the sun must give rise to the observed proton 
flux. Above that energy the CR spectrum exhibits little structure and is approximated by broken 
power laws cx E~'^: At the energy ~ 4 x 10^^ eV called the "knee", the flux of particles per 
area, time, solid angle, and energy steepens from a power law index 7 ~ 2.7 to one of index ~ 3.0. 
The bulk of the CRs up to at least that energy is believed to originate within the Milky Way 
Galaxy, typically by shock acceleration in supernova remnants. The spectrum continues with a 
further steepening to 7 ~ 3.3 at -E ~ 4 x 10^"^ eV, sometimes called the "second knee". There are 
experimental indications that the chemical composition changes from light, mostly protons, at the 
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knee to domination by iron and even heavier nuclei at the second knee 1^. This is in fact expected 
in any scenario where acceleration and propagation is due to magnetic fields whose effects only 
depend on rigidity, the ratio of charge to rest mass, Zj A. This is true as long as energy losses 
and interaction effects, which in general depend on Z and A separately, are small, as is the case 
in the Galaxy, in contrast to extra-galactic cosmic ray propagation at ultra-high energy. Above 
the so called "ankle" or "dip" at ~ 5 x 10^^ eV, the spectrum flattens again to a power law 
of index 7 ~ 2.8. This latter feature is often interpreted as a cross over from a steeper Galactic 
component, which above the ankle cannot be confined by the Galactic magnetic field, to a harder 
component of extragalactic origin. The dip at -E ~ 5 x 10^*^ eV could also be partially due to pair 
production by extra-galactic protons, especially if the extra-galactic component already starts to 
dominate below the ankle, for example, around the second-knee [H]. This latter possibility appears, 
however, less likely in light of a rather heavy composition up to the ankle suggested by several 
experiments A. In any case, an eventual cross over to an extra-galactic component is also in line 
with experimental indications for a chemical composition becoming again lighter above the ankle, 
although a significant heavy component is not excluded and the inferred chemical composition 
above ~ 10^^ eV is sensitive to the model of air shower interactions and consequently uncertain 
presently ;61. In the following we will restrict our discussion on ultra- high energy cosmic rays 
(UHECRsj above the ankle. 

Although statistically meaningful information about the UHECR energy spectrum and arrival 
direction distribution has been accumulated, no conclusive picture for the nature and distribution 
of the sources emerges naturally from the data. There is on the one hand the approximate isotropic 
arrival direction distribution [Sj which indicates that we are observing a large number of weak or 
distant sources. On the other hand, there are also indications which point more towards a small 
number of local and therefore bright sources, especially at the highest energies: First, the AGASA 
ground array claims statistically significant multi-plets of events from the same directions within a 
few degrees j71IH|, although this is controversial |H1 and has not been seen so far by the fluorescence 
experiment HiRes jlijj. The spectrum of this clustered component is oc E~^'^ and thus much 
harder than the total spectrum [7]. Second, nucleons above ~ 70EeV suffer heavy energy losses 
due to photo-pion production on the cosmic microwave background — the Greisen-Zatsepin-Kuzmin 
(GZK) effect — which limits the distance to possible sources to less than ~ 100 Mpc jl2j . For 
a uniform source distribution this would predict a "GZK cutoff" , a drop in the spectrum. However, 
the existence of this "cutoff" is not established yet from the observations and may even depend 
on the part of the sky one is looking at: The "cutoff' could be mitigated in the northern hemisphere 
where more nearby accelerators related to the local supercluster can be expected. Apart from the 
SUGAR array which was active from 1968 until 1979 in Australia, all UHECR detectors completed 
up to the present were situated in the northern hemisphere. Nevertheless the situation is unclear 
even there: Whereas a cut-off seems consistent with the few events above lO^'' eV recorded by the 
fluorescence detector HiRes ^1], it is not compatible with the 8 events above 10^" eV measured 
by the AGASA ground array ^^1- It can be remarked, however, that analysis of data based on 
a single fluorescence telescope, the so-called monocular mode in which most of the HiRes data 
were obtained, is complicated due to atmospheric conditions varying from event to event |16| . 
The solution of this problem may have to await more analysis and, in particular, the completion 
of the Pierre Auger project jl7j which will combine the two complementary detection techniques 
adopted by the aforementioned experiments and whose southern site is currently in construction 
in Argentina. 
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This currently unclear experimental situation could easily be solved if it would be possible to 
follow the UHECR trajectories backwards to their sources. However, this may be complicated 
by the possible presence of extragalactic magnetic fields, which would deflect the particles during 
their travel. Furthermore, since the GZK-energy losses are of stochastic nature, even a detailed 
knowledge of the extragalactic magnetic fields would not necessarily allow to follow a UHECR 
trajectory backwards to its source since the energy and therefor the Larmor radius of the particles 
have changed in an unknown way. Therefore it is not clear if charged particle astronomy with 
UHECRs is possible in principle or not. And even if possible, it remains unclear to which degree 
the angular resolution would be limited by magnetic deflection. This topic will be discussed in 
Sect. 3. 

The physics and astrophysics of UHECRs are also intimately linked with the emerging field 
of neutrino astronomy (for reviews see Refs. ^S|) as well as with the already established field of 
7— ray astronomy (for reviews see, e.g., Ref. JH])- Indeed, all scenarios of UHECR origin, including 
the top-down models, are severely constrained by neutrino and 7— ray observations and limits. In 
turn, this linkage has important consequences for theoretical predictions of fluxes of extragalactic 
neutrinos above about a TeV whose detection is a major goal of next-generation neutrino telescopes: 
If these neutrinos are produced as secondaries of protons accelerated in astrophysical sources and if 
these protons are not absorbed in the sources, but rather contribute to the UHECR flux observed, 
then the energy content in the neutrino flux can not be higher than the one in UHECRs, leading to 
the so called Waxman-Bahcall bound for transparent sources with soft acceleration spectra |2n|l21j. 
If one of these assumptions does not apply, such as for acceleration sources with injection spectra 
harder than E~'^ and/or opaque to nucleons, or in the top-down scenarios where X particle decays 
produce much fewer nucleons than 7— rays and neutrinos, the Waxman-Bahcall bound does not 
apply, but the neutrino flux is still constrained by the observed diffuse 7— ray flux in the GeV 
range. This will be discussed in the following section. 

2 Severe Constraints on Scenarios producing more photons than 
hadrons 

Electromagnetic (EM) energy injected above the threshold for pair production on the cosmic mi- 
crowave background (CMB) at ~ 10^^/(1 + z)eV at redshift z (to a lesser extent also on the 
infrared/optical background, with lower threshold) leads to an EM cascade, an interplay between 
pair production followed by inverse Compton scattering of the produced electrons. This cascade 
continues until the photons fall below the pair production threshold at which point the universe 
becomes transparent for them. In todays universe this happens within just a few Mpc for injection 
up to the highest energies above 10^*^ eV. All EM energy injected above ~ 10^^ eV and at distances 
beyond a few Mpc today is therefore recycled to lower energies where it gives rise to a characteristic 
cascade spectrum oc E~'^'^ down to fractions of a GeV 22, . The universe thus acts as a calorimeter 
where the total EM energy injected above ~ 10^^/(1 + z) eV is measured as a diffuse isotropic 
7— ray flux in the GeV regime. This diffuse flux is not very sensitive to the somewhat uncertain 
infrared/optical background '231. Any observed diffuse 7— ray background acts as an upper limit 
on the total EM injection. Since in any scenario involving pion production the EM energy fluence 
is comparable to the neutrino energy fluence, the constraint on EM energy injection also constrains 
allowed neutrino fluxes. 
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This diffuse extragalactic GeV 7— ray background can be extracted from the total 7— ray flux 
measured by EGRET by subtracting the Galactic contribution. Since publication of the original 
EGRET limit in 1995 models for this high latitude Galactic 7— ray foreground were improved 
significantly. This allowed the authors of Ref. [^5] to reanalyze limits on the diffuse extragalactic 
background in the region 30 MeV-10 GeV and to lower it by a factor 1.5-1.8 in the region around 
1 GeV. There are even lower estimates of the extragalactic diffuse 7— ray flux j^HI- In this article, 
however, we will use the more conservative limits from Ref. |25| . 

The energy in the extra-galactic 7— ray background estimated in Ref. ^25^ is slightly more than 
one hundred times the energy in UHECR above the GZK cutoff. The range of such trans-GZK 
cosmic rays is about ~ 30Mpc, roughly one hundredth the Hubble radius, and only sources within 
that GZK range contribute to the trans-GZK cosmic rays. Therefore, any mechanism involving 
sources distributed roughly uniformly on scales of the GZK energy loss length ~ SOMpc and 
producing a comparable amount of energy in trans-GZK cosmic rays and photons above the pair 
production threshold can potentially explain this energy flux ratio. The details depend on the exact 
redshift dependence of source activity and other parameters and in general have to be verified by 
numerically solving the relevant transport equations, see, e.g., Ref. |2Z1- Such mechanisms include 
shock acceleration in powerful objects such as active galactic nuclei j28j . 

On the other hand, any mechanism producing considerably more energy in the EM channel 
above the pair production threshold than in trans-GZK cosmic rays tend to predict a ratio of the 
diffuse GeV 7— ray flux to the trans-GZK cosmic ray flux too high to explain both fluxes at the 
same time. As a consequence, if normalized at or below the observational GeV 7— ray background, 
such scenarios tend to explain at most a fraction of the observed trans-GZK cosmic ray flux. Such 
scenarios include particle physics mechanisms involving pion production by quark fragmentation, 
e.g. extra-galactic top-down mechanisms where UHECRs are produced by fragmenting quarks re- 
sulting from decay of superheavy relics [20] • Most of these quarks would fragment into pions rather 
than nucleons such that more 7— rays (and neutrinos) than cosmic rays are produced. Overpro- 
duction of GeV 7— rays can be avoided by assuming the sources in an extended Galactic halo with 
a high > 10'^ overdensity compared to the average cosmological source density, which would also 
avoid the GZK cutoff [HOI- These scenarios, however, start to be constrained by the anisotropy they 
predict because of the asymmetric position of the Sun in the Galactic halo for which there are no 
indications in present data [HJ. Scenarios based on quark fragmentation also become problematic 
in view of a possible heavy nucleus component and of upper limits on the photon fraction of the 
UHECR flux 0. 

As a specific example for scenarios involving quark fragmentation, we consider here the case 
of decaying Z-bosons. In this "Z-burst mechanism" Z-bosons are produced by UHE neutrinos 
interacting with the relic neutrino background (22j- If the relic neutrinos have a mass rriiy, Z-bosons 
can be resonantly produced by UHE neutrinos of energy E^, ~ M^; / {2my) ~ 4.2 x 10^^ eV {eV/rriu)- 
The required neutrino beams could be produced as secondaries of protons accelerated in high- 
redshift sources. The fluxes predicted in these scenarios have recently been discussed in detail, 
for example, in Refs. |33[ I27j . In Fig. ^ we show an optimistic example taken from Ref. It 
is assumed that the relic neutrino background has no significant local overdensity. Furthermore, 
the sources are assumed to not emit any 7— rays, otherwise the Z-burst model with acceleration 
sources overproduces the diffuse GeV 7— ray background 1131 • We note that no known astrophysical 
accelerator exists that meets the requirements of the Z-burst model |43| I44j . 

However, a combination of new constraints allows to rule out that the Z-burst mechanism 



4 



explains a dominant fraction of the observed UHECR flux, even for pure neutrino emitting sources: 
A combination of cosmological data including the WMAP experiment limit the sum of the masses of 
active neutrinos to < 1 eV (45. . Solar and atmospheric neutrino oscillations indicate that individual 
neutrino masses are nearly degenerate on this scale ,46,, and thus the neutrino mass per flavor 
must satisfy < 0.33 eV. However, for such masses phase space constraints limit the possible 
over-density of neutrinos in our Local Group of galaxies to < 10 on a length scale of ~ 1 Mpc 47 . 
Since this is considerably smaller than the relevant UHECR loss lengths, neutrino clustering will 
not significantly reduce the necessary UHE neutrino flux compared to the case of no clustering. For 
the maximal possible value of the neutrino mass rrii, ~ 0.33 eV, the neutrino flux required for the 
Z-burst mechanism to explain the UHECR flux is only in marginal conflict with the FORTE upper 
limit , and factor 2 higher than the new GLUE limit , as shown in Fig. ^ For all other 
cases the conflict with both the GLUE and FORTE limits is considerably more severe. Also note 
that this argument does not depend on the shape of the low energy tail of the primary neutrino 
spectrum which could thus be even mono-energetic, as could occur in exclusive tree level decays of 
superheavy particles into neutrinos |1H1- However, in addition this possibility has been ruled out 
by overproduction of GeV 7— rays due to loop effects in these particle decays 1^ . 

The possibility that the observed UHECR flux is explained by the Z burst scenario involving 
normal astrophysical sources which produce both neutrinos and photons by pion production is 
already ruled out by the former EGRET limit: In this case the GeV 7— ray flux level would have 
roughly the height of the peak of the neutrino flux multiplied with the squared energy in Fig. ^ 
thus a factor ~ 100 higher than the EGRET level. 

Any further reduction in the estimated contribution of the true diffuse extra-galactic 7— ray 
background to the observed flux, therefore, leads to more severe constraints on the total EM 
injection. For example, future 7— ray detectors such as GLAST I50_ will test whether the diffuse 
extragalactic GeV 7— ray background is truly diffuse or partly consists of discrete sources that could 
not be resolved by EGRET. Astrophysical discrete contributions such as from intergalactic shocks 
are in fact expected This could further improve the cascade limit to the point where even 

acceleration scenarios may become seriously constrained. 

3 Cosmic Magnetic Fields and Their Influence on Ultra-High En- 
ergy Cosmic Ray Propagation 

Cosmic magnetic fields are inextricably linked with cosmic rays in several respects. First, they play a 
central role in Fermi shock acceleration. Second, large scale extra-galactic magnetic fields (EGMF) 
can cause significant deflection of charged cosmic rays during propagation and thus obviously 
complicate the relation between observed UHECR distributions and their sources. 

Magnetic fields are omnipresent in the Universe, but their true origin is still unclear [^2]. Mag- 
netic fields in galaxies are observed with typical strengths of a few micro Gauss, but there are also 
some indications for fields correlated with larger structures such as galaxy clusters |53| . Magnetic 
fields as strong as ~ IfiG in sheets and filaments of the large scale galaxy distribution, such as in 
our Local Supercluster, are compatible with existing upper limits on Faraday rotation |531 1541 . 
It is also possible that fossil cocoons of former radio galaxies, so called radio ghosts, contribute 
significantly to the isotropization of UHECR arrival directions j56j . 

To get an impression of typical deflection angles one can characterize the EGMF by its r.m.s. 
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strength B and a coherence length Ic- If we neglect energy loss processes for the moment, then the 
r.m.s. deflection angle over a distance r > Ic in such a field is 9{E,r) ~ {2rlc/9)^^'^ /tl |57j . where 
the Larmor radius of a particle of charge Ze and energy is r/^ ~ i?/ {ZeB). In numbers this reads 

for r > Ic- This expression makes it immediately obvious that fields of fractions of micro Gauss 
lead to strong deflection even at the highest energies. This goes along with a time delay t{E, r) ~ 
re{E, df/A ~ 1.5 X 10^ Z'^{E/l{)^^ eV)-2(r/10 Mpc)2(/ jMpc)(5/10-9 G)^ yr which can be millions 
of years. A source visible in UHECRs today could therefore be optically invisible since many models 
involving, for example, active galaxies as UHECR accelerators, predict variability on shorter time 
scales. 

Quite a few simulations of the effect of extragalactic magnetic fields (EGMF) on UHECRs 
exist in the literature, but usually idealizing assumptions concerning properties and distributions 
of sources or EGMF or both are made: In Refs. |58| I59| I6fll 1611 162) sources and EGMF follow 
a pancake profile mimicking the local supergalactic plane. In other studies EGMF have been 
approximated in a number of fashions: as negligible [63| I64| . as stochastic with uniform statistical 
properties |65l 1661 16 7j , or as organized in spatial cells with a given coherence length and a strength 
depending as a power law on the local density jBEl- Only recently attempts have been made 
to simulate UHECR propagation in a realistically structured universe j69| I71j . For now, these 
simulations are limited to nucleons. 

In Ref. jnn] the magnetized extragalactic environment used for UHECR propagation is produced 
by a simulation of the large scale structure of the Universe. The simulation was carried out within 
a computational box of 50 Mpc length on a side, with normalized Hubble constant h = i?o/(100 
km s^^ Mpc~^) = 0.67, and using a comoving grid of 512^ zones and 256'^ dark matter particles. 
The EGMF was initialized to zero at simulation start and subsequently its seeds were generated at 
cosmic shocks through the Biermann battery mechanism jTU]. Since cosmic shocks form primarily 
around collapsing structures including filaments, the above approach avoids generating EGMF in 
cosmic voids. 

In Ref. j2J constrained simulations of the local large scale structure were performed and the 
magnetic smoothed particle hydrodynamics technique was used to follow EGMF evolution. The 
EGMF was seeded by a uniform seed field of maximal strength compatible with observed rotation 
measures in galaxy clusters. 

The questions considered in these two works were somewhat different, however. In Ref. [7J 
deflections of UHECR above 4 x 10^^ eV were computed as a function of the direction to their source 
which were assumed to be at cosmological distances. This made sense, because (i) the constrained 
simulations gives a viable model of our local cosmic neighborhood within about 100 Mpc, at least 
on scales beyond a few Mpc and (ii) the deflections typically were found to be smaller than a few 
degrees. Concrete source distributions were not considered. 

In contrast, Ref. fBP_ was not concerned with concrete sky distributions or deflection maps 
because the simulation was unconstrained and thus only gave a typical large scale structure model 
and not our concrete local neighborhood. Instead, the question was asked which observer positions 
and source distributions and characteristics lead to UHECR distributions whose spherical multi- 
poles for / < 10 and auto-correlation at angles 9 < 20° are consistent with observations. As a result 
it was found that (i) the observed large scale UHECR isotropy requires the neighborhood within 
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a few Mpc of the observer is characterized by weak magnetic fields below 0.1 fiG, and (ii) once 
that choice is made, current data do not strongly discriminate between uniform and structured 
source distributions and between negligible and considerable deflection. Nevertheless, current data 
moderately favor a scenario in which (iii) UHECR sources have a density ~ 10~^ Mpc^^ and 
follow the matter distribution and (iv) magnetic fields are relatively pervasive within the large 
scale structure, including filaments, and with a strength of order of a /iG in galaxy clusters. A 
two-dimensional cut through the baryonic density and EGMF environment of the observer in a 
typical such scenario is shown in Fig. [2 

It was also studied in Ref. [22] how future data of considerably increased statistics can be 
used to learn more about EGMF and source characteristics. In particular, low auto-correlations 
at degree scales imply magnetized sources quite independent of other source characteristics such 
as their density. The latter can only be estimated from the auto-correlations halfway reliably if 
magnetic fields have negligible impact on propagation. This is because if sources are immersed 
in considerable magnetic fields, their images are smeared out, which also smears out the auto- 
correlation function over several degrees. For a sufficiently high source density, individual images 
can thus overlap and sensitivity to source density is consequently lost. The statistics expected from 
next generation experiments such as Pierre Auger ^7] and EUSO [721 should be sufficient to test 
source magnetization by the auto-correlation function |69j . 

Interestingly, however, there is a considerable quantifiable difference in the typical deflection 
angles predicted by the two EGMF scenarios in Refs. |691 171j that can not be compensated by 
specific source distributions: Even for homogeneous source distributions, the average deflection 
angle for UHECRs above 4 x 10^^ eV obtained in Ref. [22] is much larger than in Ref. [JJ, as 
can be seen in Fig. |31 In fact, even if the magnetic field strength is reduced by a factor 10 in the 
simulations of Ref. I69j . the average defiection angle above 4 x 10^^ eV is still ~ 30°, only a factor 
~ 2.2 smaller. This non-linear behavior of defiection with field normalization is mostly due to the 
strongly non-homogeneous character of the EGMF. 

Most of these differences are probably due the different numerical models for the magnetic fields. 
Although Ref. [7J start with uniform seed fields, whereas in Ref. j69j seed fields are injected at 
shocks, by itself, this difference should not infiuence the resulting EGMF very much at late times, 
at least inside galaxy clusters [31] . It should be noted, however, that in the filaments, where the 
gas motions are more uniform, the simulated magnetic fields may depend to a certain extent on the 
initial seed fields although that is not trivial to quantify in general terms. In addition, numerical 
resolution may play an important role because it affects the amplification and the topological 
structure of the magnetic fields, both of which are important for the normalization procedure, see 
below. The resolution in Ref. 'M' is constant and much better in filaments and voids but worse 
in the core of galaxy clusters than the (variable) resolution in Ref. |71j . If in both simulations the 
magnetic fields are normalized to (or reproduce) the same "observed" values in the core of rich 
clusters then obviously their values in the filaments will be very different for the reasons outlined 
above. This may partly explain why the contribution of filaments to UHECR deflection is more 
important in Ref. although a more detailed analysis and comparison are required to settle the 
issue. In any case, the magnetic fields obtained in Ref. [5^] seem to be quite extended, as can be 
seen in Fig. ^ About 10% of the volume is filled with fields stronger than 10 nano Gauss, and a 
fraction of 10~^ is filled by fields above a micro Gauss. The different amounts of defiection obtained 
in the simulations of Refs. |69|l71j show that the distribution of EGMF and their effects on UHECR 
propagation are currently rather uncertain. 
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Finally we note that these studies should be extended to include heavy nuclei [73] since there 
are indications that a fraction as large as 80% of iron nuclei may exist above 10^^ eV [H]. As a 
consequence, even in the EGMF scenario of Ref. [7J deflections could be considerable and may 
not allow particle astronomy along many lines of sight: The distribution of deflection angles in 
Ref. j71j shows that deflections of protons above 4 x 10^^ eV of > 1° cover a considerable fraction 
of the sky. Suppression of deflection along typical lines of sight by small filling factors of deflectors 
is thus unimportant in this case. The deflection angle of any nucleus at a given energy passing 
through such areas will therefore be roughly proportional to its charge as long as energy loss lengths 
are larger than a few tens of Mpc [3]. Deflection angles of ~ 20° at ~ 4 x 10^^ eV should thus 
be the rule for iron nuclei. In contrast to the contribution of our Galaxy to deflection which can 
be of comparable size but may be corrected for within sufficiently detailed models of the galactic 
field, the extra-galactic contribution would be stochastic. Statistical methods are therefore likely to 
be necessary to learn about UHECR source distributions and characteristics. In addition, should 
a substantial heavy composition be experimentally confirmed up to the highest energies, some 
sources would have to be surprisingly nearby, within a few Mpc, otherwise only low mass spallation 
products would survive propagation [7H] . 

The clustered component of the UHECR spectrum may play a key role in this context and 
may be caused by discrete sources in directions with small deflection. Since, apart from energy 
losses, cosmic rays of same rigidity Zj A are deflected similarly by cosmic magnetic fields, one may 
expect that the composition of the clustered component may become heavier with increasing energy. 
Indeed, in Ref. [TBj it was speculated that the AGASA clusters may be consistent with consecutive 
He, Be-Mg, and Fe bumps. 

4 Conclusions 

We have reviewed two current issues in theoretical ultra-high energy cosmic ray research. 

The first one concerns constraints on scenarios attempting to explain highest energy cosmic rays 
by extra-galactic sources producing not only cosmic rays but also photons: Improved data analysis 
and, in the future, improved data for example from GLAST can considerably reduce estimates 
of the true extra-galactic GeV 7— ray background which acts as a calorimeter of electromagnetic 
energy injected above ~ 10^^/ (1 + 2;) eV. Already current estimates imply that scenarios producing 
considerably more photons than hadrons, such as extra-galactic top-down scenarios and the Z-burst 
mechanism, can not explain all of the highest energy cosmic ray flux. A further reduced diffuse 
GeV 7— ray background will start to constrain even normal acceleration scenarios. 

As for the second issue we pointed out that the influence of large scale cosmic magnetic fields 
on ultra-high energy cosmic ray propagation is currently hard to quantify and may not allow to do 
"particle astronomy" along most lines of sight, especially if a significant heavy nucleus component 
is present above 10^^ eV. In this case extensive Monte Carlo simulations including nuclei and based 
on constrained large scale structure simulations will be necessary to fully exploit data from future 
instruments such as the Pierre Auger Wj\ and EUSO projects |72j . 
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E[eV] 

Figure 1: Flux predictions for a Z-burst model averaged over flavors and characterized by a neutrino 
injection flux per comoving volume oc up to 3 x W^^ eV and for redshifts between and 3. The 
sources are assumed to be exclusive neutrino emitters. All neutrino masses were assumed equal with 
rui, = 0.33 eV and we again assumed maximal mixing between all flavors. Also shown are predicted 
and observed cosmic ray and 7— ray fluxes, the atmospheric neutrino flux |34j . as well as existing 
upper hmits on the diffuse neutrino fluxes from MACRO AMANDA II jSHl, BAIKAL |37j, 
AGASA 38,, the Fly's Eye ^ and RICE HQ] experiments, and the limits obtained with the 
Goldstone radio telescope (GLUE) [31] and the FORTE satellite [32], as indicated. The cosmic ray 
data are from the AGASA JHj and HiRes Jl] experiments, and the new EGRET estimate of the 
extra-galactic diffuse 7— ray flux is shown to the left. From Ref. j27j . 
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Figure 2: Log-scale two-dimensional cut through magnetic field total strength in Gauss (color scale 
in Gauss, upper panel) and baryon density in units of average baryon density (color scale, lower 
panel), for a scenario in good agreement with UHECR data studied in Ref. j69j . The observer is in 
the center of the figures and is marked by a star. The EGMF strength at the observer is ~ 10~^^ G. 
Note that both panels correspond to the same cuts through the full large scale simulation box used. 
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Figure 3: The cumulative distribution of minimal UHECR deflection angles a with respect to the 
line of sight to the sources. This is for a scenario from Ref. j69j in good agreement with UHECR 
data, where the sources follow the baryon density and have average density Ug = 2.4 x 10~^ Mpc~'^, 
and the EGMF included in the large scale structure simulation reaches several micro Gauss in 
the most prominent galaxy cluster. Shown are the average (middle, histogram) and 1-a variations 
(upper and lower curves) above 4 x 10^^ eV, over 24 realizations varying in the positions and 
luminosities Qi of individual sources, the latter assumed to be distributed as dus/dQi oc ' with 
1 < Qi < 100 in arbitrary units. Also given on top of the figure are average and variances of the 
distributions. 
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Figure 4: The cumulative filling factors for EGMF strength in the simulations used in Ref. 
above (decreasing curve) and below (increasing curve) a given threshold, as a function of that 
threshold. 
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